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Poincare gauge theory (PGT) is an alternative gravity theory, which attempts to bring the gravity 
into the gauge-theoretic frame, where the Lagrangian is quadratic in torsion and curvature. Recently, 
the cosmological models with torsion based on this theory have drawn many attentions, which try 
to explain the cosmic acceleration in a new way. Among these PGT cosmological models, the one 
with only even parity dynamical modes - SNY model, for its realistic meaning, is very attractive. 
In this paper, we first analyze the past-time cosmic evolution of SNY model analytically. And based 
l" ■ on these results we fit this model to the most comprehensive SNela data (Union 2) and thus find the 

I ' best-fit values of model parameters and initial conditions, whose related value is consistent with 

, the one from ACMD at the la level. Also by the estimate, we provide certain constraints on 

. these parameters. Using these best-fit values for the Union 2 SNela dataset, we are able to predict 

' the evolution of our real universe over the late time. From this prediction, we know the fate of our 

^ [ universe that it would expand forever, slowly asymptotically to a halt, which is in accordance with 

, ^ ■ the earlier works. 

o ■ 

U: 

About 12 years ago, the high redshift SNela observations, suggesting that our universe is not only expanding but 
5—1 ' also accelerating, have changed our basic understanding of the universe, which led us to the new era of cosmology. In 
order to explain this weird phenomena, cosmologists introduced the cosmological constant back again, and constructed 
' a new standard model also known as the concordance model, which is intended to satisfy all the main observations 
T— I such as type la supernovae (SNela), cosmic microwave background radiation (CMBR) and large scale structure (LSS). 
J> Although the cosmological constant accounts for almost 74% of the whole energy density in this concordance model, 
, the value is still too small to be explained by any current fundamental theories. Lacking the underlying theoretical 
9^ ' foundations, the particular value of cosmological constant is just selected phenomcnologically, which means the model 
, is highly sensitive to the value of model parameter, resulting in the so-called the fine-tuning problem. This problem 
• ■ is considered as the biggest issue for almost all cosmological models. In order to alleviate this troublesome problem, 
|, various dynamical dark energy theories have been proposed and developed these years, such as quintessence [l|, Q and 
^ I ■ phantom ^^i^, in which the energy composition depends on time. But these exotic fields are still phenomenological, 
I I lacking theoretical foundations. Besides adding some unknown fields, there is another kind of theories known as 
• • . modified gravity, which use alternativegravity theory instead of Einstein theory, such as f{R) theory [1, Q, MOND 
' cosmology 7], Poincare gauge theory [8|,|1], and de Sitter gauge theory [iQ. Among these theories, Poincare gauge 
K> theory (PGT) has a solid theoretical motivation, which tries to bring the gravity into the gauge-theoretic framework. 
^ The attempt to treat gravitation as a gauge interaction could date back to the 1950s. Utiyama first presented 

^ \ a groundwork for a gauge theory of gravitation [llj. Later Kibble and Sciama [l2| inherited Utiyama's work and 
developed a gauge theory of gravitation commonly called Einstein-Cartan-Sciama-Kibble (ECSK) theory, where the 
symmetry group is Poincare group. In ECSK theory, the Lagrangian has only one linear curvature term, thus the 
equation of torsion is algebraic, and torsion is non-dynamical. Some follow-up cosmological models with torsion, based 
on ECSK model, were investigated, which try to avoid the cosmological singularity 13], and from these cosmological 
discussions, they found the torsion was imagined as playing role only at high densities in the early universe. In order 
to enable torsion to propagate, one has to introduce quadratic term of torsion and curvature, which is spoken of 
Poincare gauge theory, proposed by Hehl in 1980 (l^ . Poincare group is a semiproduct of translation group and 
Lorentz rotation group, which is the global symmetry group of a space-time in absence of gravity. If one localize 
this global symmetry as a gauge symmetry, the gravity would emerge spontaneously. Note that in this transition 
some independent compensating fields have to be introduced, the orthonormal coframe and the metric-compatible 
connection, which correspond to the translation and local rotation potentials, respectively. And their corresponding 
field strength are torsion and curvature. The early works of PGT are intended to unify gravity and the other three 
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interactions in the framework of gauge field theory. However now it might be the solution to the cosmic acceleration 
problem, in which the dynamical torsion could play the role as dark energy to drive our universe to accelerate [Tsl - fisj . 

In PGT, the propaga ting torsion field can be identified as six possible dynamic modes, carrying spins and parity: 
2^, 1^,0^, respectively |l9l [20|. However, expect for two "scalar modes", 0^, other dynamical modes are not phys- 
ically acceptable, which are ruled out by certain theoretical constraints (e.g., "no-ghosts" and "no-tachyons"). The 
pseudoscalar mode with the odd parity, 0^, is reflected in the axial vector torsion, which is driven by the intrinsic 
spin of elementary fermions. It is generally thought that the axial torsion is large and has notable effects only in the 
very early universe, where the spin density is large. Therefore, this part probably have insignificant contribution to 
the current evolution of our universe. Whereas the scalar mode with the even parity, 0+, which is associated with the 
vector torsion, does not interact in any direct way with any known of matter. Consequently, it is rational to imagine 
it has an significant magnitude but yet not been noticed. In 2008, Shie, Nester and Yo (SNY) proposed a cosmological 
model based on PGT, where the dynamical scalar torsion 0+ accounts for the cosmic acceleration lip. Ilal . Later, 
we investigated the dynamical properties of SNY model and applied the statefinder diagnostics to it[2l|, [23, and did 
some analytical analysis on the late-time evolution in [2^. In 2010, Baekler, Hehl and Nester generalized the SNY 
model to the BHN model which includes the 0~ mode and even the cross parity couplings and gave a comprehensive 
picture of PGT cosmology And some follow-up works suggest the cosmic acceleration is still mainly due to the 
even mode. [13, [3 

In this paper, we study the evolution of our universe based on SNY model. In SNY model, the cosmic evolution 
is described in terms of three variables, {H,^,R), by three first-order dynamical equations. We obtain three exact 
solutions of scale factor with certain particular constant affine scalar curvature, and show that these solutions are 
not physically real, which requires us to study the non-constant curvature case. Then we extend the investigation 
to the non-constant affine scalar curvature, and find the analytical solutions for the past-time evolution, which could 
be compared with the observational data to constrain the values of parameters and determine the goodness of this 
model. Then we fit these theoretical results to the Union 2 supernovae dataset, and find the best-fit values of the model 
parameters and initial conditions, by means of the estimate. Also, we plot the contours of certain confidence levels, 
which shows some constraints on the parameters. Finally, we study the future evolution. Using these constraints of 
parameters, we predict that the real universe would expand forever, slowly asymptotically to a halt, which is consistent 
with some earlier works. p3. [23j 



II. SNY MODEL AND ITS SOME EXACT SOLUTIONS 

In order to localize the global Poincare symmetry, one has to introduce two compensating fields: coframe field 1?" 
and metric-compatible connection field F"^ — F^'^'^'da:* ^. And their associated field strengths are the torsion and 
curvature 2-forms 

= ir^^ 1?" A i?^ = d'^^' + Fj; A i?'' (1) 
i?^'^ = ii?;^^ I?" Ai9^ ^dF^'^-FF^AF"" (2) 

which satisfy the respective Bianchi identities: 

DT" = mi A r, DR1 = (3) 

The Lagrangian denstiy in PGT takes the standard quadratic Yang-Mills form, qualitively, 

F] ~ curvature -I- (torsion)^ + (curvature)^. (4) 

Some early works on PGT concluded that in the weak field approximation, the torsion field could be identified as 
six irreducible parts with certain dynamical modes, which propagate with 2^, 1*, 0*. Later some investigations by 
linearized theory and Hamiltonian analysis concluded that only the two "scalar modes" are physically acceptable, 
carrying 0^. Furthermore, because the 0~ is driven by spin density, this mode would not have a significant effects, 
expect in the very early universe. And some numerical results demonstrate that the acceleration is mainly due to the 



^ the Greek indices a, /9, 7... denote the 4d orthonormal indices, whereas the Latin indices denote 4d coordinate (holonomic) indices and 
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0+. So here we only consider the simple 0+ case, i.e., SNY model. In this case, the gravitational lagrangian density 
is simplified to this specific form. 



-aoR + Y,<^nT^ + —R 



(5) 



where is the algebraically irreducible parts of torsion, R is the scalar curvature, which is the non- vanishing irreducible 
parts of curvature [24] . Note that ao and an are dimensionless parameters, whereas b have the same dimension with 

R-\ 

Since current observations favor a homogeneous, isotropic and spatially flat universe, it is reasonable to work on 
the Friedman-Lemaitre-Robertsen- Walker (FLRW) metric, where the isotropic orthonormal coframe takes the form: 



= dt, ^ a{t)dx'; 

and the only non- vanishing connection 1-form coefficients are of the form: 

Consequently, the nonvanishing torsion tensor components take the form 



(6) 
(7) 

(8) 



where $ represents torsion. By variation w.r.t. the coframe and connection, one could find the cosmological equations 
of SNY model. 
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where the H = a/ a is the Hubble parameter, /i = 02 — oq, which represents the mass of 0+ mode, and the energy 
density of matter component is 
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(12) 



Note that the universe here is assumed to be the dust universe, for the effect of radiation is almost negligible in the 
late-time evolution. The Newtonian limit requires ao = —1. 

From Eq. (jll[) , it is easy to find the scalar afhne curvature remains a constant R — ~6fi/b forever as long as its 
initial data has this special value. In this case, Eq. ([9]) can be rewritten as 



2ad + + --^a^ = 0, 
2 a20 



(13) 



where the scale factor a{t) is decoupled to $ and R fields. Thus, we could obtain some simple exact solution of a{t). 
The positivity of the kinetic energy requires 02 > and 6 > [2^ , so we have the solution 



a(t) = ao 



cos 


^ (^ — ^o) — arctan 






COS 


arctan 







(14) 



where ^ = iJ,/\/2a2b and Hq = H{to). However, such a choice confiicts with the assumption of energy positivity in 
the R — —6fi/b case. 

If we audaciously relax the parameter requirement for positive kinetic energy, i.e., 02 < —1 and fi < 0, this phantom 
scenario will turn out to be interesting. Now we have the solution 



a{t) — ao exp 



(g-2)(t-M + Aln(^'^°) + (^ + ^°)---P[^^(^-^°)] 
3^ 2^ 



(15) 



4 



where ^ = fx/ V— 2a2&. From this expression, it is obvious that the late-time behavior would be analogous with the 
exponential expansion of the de Sitter universe. In other words, dark energy can be mimicked in such case. Using the 
dynamical analysis, we have pointed out that there is a late-time de Sitter attractor [111, [2^. Note that the solution 
([T5|) is just corresponding to the de Sitter attractor. 
Especially, as 02 = —1, we have a solution 

From this expression, we could find that the late-time evolution is similar to the t'^^^ expansion of the matter dominant 
universe. 

In this constant affine scalar curvature case, the three exact solutions we find are all not physically acceptable: the 
first conflicts with the assumption of energy positivity; while the second and third ones breaks the positivity of kinetic 
energy. For this reason, these solutions with constant scalar curvature cannot describe our real universe, and thus we 
have to study the case of non-constant scalar curvature. 

Unfortunately, for such a complex nonlinear system, Eqs. (|9))-([TT|). the whole evolution of our universe is difficult to 
obtain, so we are forced to divide it into two seperate parts, the past and the future, and discuss them respectively. 
Among these two parts, it is the past one that could be compared with the existing observational results and impose 
constraints on the values of model parameters and initial conditions. Based on these constraints, we could obtain the 
future evolution and the fate of our real universe. Therefore, we study the past-time evolution first. 



III. THE ANALYTICAL APPROACH FOR THE a{t) < a{to) 

In Ref.[23|, we analyzed the evolution when a{t) > a(to), i.e. the future evolution. In this section, we investigate 
the evolution when a{t) < a(to), i-e. the evolution of universe at t < Iq. 
For such a dust universe, the continuity equation is 

P = -SHp, (17) 

which could also be derived directly by the dynamical equations Eqs. ([^- (ITT]) and the constraint equation Eq. (|I2I) . 
This equation is easy to solve, 

P=^, (18) 

where po is the current matter density of our universe. 
If we rescale the variables and parameters as 

t^t/lo; H^IqH; R ^ l^R; 

$^^0*; P^llup; b^b/ll (19) 

where Iq = c/Hq is the Hubble radius, these motion equations Egs.ipil- dTT]) would be dimensionless. 

By the substitution Eq. p8)) back to Eas.([9))- (ITT]) and rescaling transformations, one could obtain the new equation 
set, which is more convenient to solve, 

aHH' = JLR^f^L^2H^ (20) 
602 602 a'^ 



aH<^' = -^i-_3H$+i$2 (21) 

2a2 2a2 3 ^ ' 



aHR' - -^(i?+^)ci>, (22) 
3 

where the prime means the derivation with respect to scale factor a. The current scale factor is generally supposed 
to be unity, so it is reasonable to assume the ansatz is as follows. 



H{a) = ho + Y.K{a-lT, (23) 

00 

$(a) = ^o + ^'/5«(a-l)", (24) 

n=l 

00 

R{a) = ro + ^r„(a-l)", (25) 



n=l 
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which has a good convergence when a{t) < a(io).^ 

Substitute this ansatz back to the Eg s. ([201) - ([2^ . we find the recursion relation of cocfiicicnts, 



when n > 2, 



with 



hi 

ri -- 



J_ / MH) _ _Po_ 
ho \6a2 6a2 

J_ [ II. _ Pi. 
ho \2a2 2a2 



2ht 



Shoipo 



(26) 
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(28) 
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i=0 
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1 
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(29) 
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4a2 — ' ^ — ' 



n-l 



i=0 



1=0 
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1=0 
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1 
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1=0 
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(30) 
(31) 

(32) 



where /iq, and tq are the initial conditions of iJ(l), $(1) and -R(l), i.e. the present time. Note that for the rescaling 
transformation Eg. dTO)) . the initial value for iJ is supposed to be unity, and therefore ft-o = 1- We here obtained the 
analytical solution of the past evolution of SNY model, which can be tested by the observational results, realistically. 



IV. THE CONSTRAINTS OF PARAMETERS VIA SNEIA 



The most common approach to test a cosmological model is the supernovae fitting. In this section we attempt to 
fit the model parameters and the initial values to the current type la supernovae data. 

The supernovae data we use here is the Union 2 dataset (N=557), which is the most comprehensive one up to date, 
combining the former SNela dataset in a homogeneous manner. It consists of distance modulus /itobs, which equals 
to the difference between apparent magnitude rtii and the absolute magnitude M^, redshifts Zi of supernovae, and 
the covariance matrix Csn represents the statistical and systematic errors. By comparing the theoretical distance 
modulus [ith-, derived from the cosmological model, with the observational data, we could obtain the constraints on 
the model parameters and initial value as well as the goodness of the model. 

As stated above, the SNY model predicts a specific form of the Hubble constant H(a;a2,b,(po,ro) as a function 
of scale factor, which is explicitly expressed in Eqs. (P^ . (l251) and 1^^. For the relation between scale factor a and 
redshift z, it is easy to rewrite H(a) in terms of redshift. 



H{z;a2,b,(po,ro) = ^(-l)"/i„ 



ri=0 



(33) 



^ Here, a{t) is set to unity. 
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which is more convenient to compare to the supernovae data. 

The theoretical distance modulus is related to the luminosity distance by 



tJ-th{zi) = 5 log. 



10 



Mpc 



25 



Slogio DL{zi) - 51ogi 



V Mpc 

= 51ogioi?L(20-51ogio/i + 42.38 



25 



(34) 



where the Dl{z) is the dimensionless 'Hubble-constant free' luminosity distance defined by Dl{z) = HodL{z)/c. 

For the spatially flat model we consider here, the 'Hubble-constant free' luminosity distance could be expressed in 
terms of Hubble parameter H{z; 02, 6, 4>o,ro), Eq. 



Dl{z) = (1 + 2) / dz' 



1 



H{z']a2,b,(l)o,ro)' 



(35) 



Due to the normal distribution of errors, the method could be used as the maximum likelihood estimator to 
compare the theoretical models and the observational data, which would show us the best-fit parameters (02 , 6, tpo 
, ro) and the goodness of model. The here for the SNela data is 



N 



hi 

N 

= X! i^^obs(zi) - 51ogio DL{zi;0) - ^lo] {Cg^)ij[fj.obs{zj) - 5 log^g L»l(zj; 6*) - /iq]- 



(36) 



where /iq — — 51og]^o ^ + 42.38 and 9 denotes the model parameters and initial values. The parameter fiQ is a nuisance 
parameter, whose contribution we are not interested in. By marginalization over this parameter ^q, we obtain 



c_ 

271 I ' 



where 



N 



M^) = i^^obs{z^;e) -5\og^QDL{z^;e)] {Cs^).,j[fiobs{zj;e) -5\og,oDL{z.r,e)], 

N 

B{e) = Y,{Csl)^,\^J^obs{z,■6)-h\0g^^DL{z^M. 

3 

N 



(37) 

(38) 
(39) 
(40) 



By minimizing the x^, we find the best- fit values of parameters and initial conditions (02, 6, (po, ?'o) of SNY model ^, as 
shown in Tabffl The minimal here is 535.284 , whereas the value for ACDM is 536.634, with VLm = 0.27, I^a = 0.73, 



a2 b (fo ro 




1.336 0.992 0.584 5.839 


535.284 



TABLE I: The best-fit initial data and parameters 
which implies that the torsion cosmology is consistent with ACMD at the la level. Substitute these best-fit values 



^ The full numerical calculation are performed using Matlab here. 



7 



1.5 



1.4- 




0.8- 



0.7 1 • • • • 1 

1 1.5 2 2.5 

FIG. 1: The 68.3%, 95.4% and 99.7% confidence contours with respect to model parameters 02, b, using the Union 2 
dataset. Here we assume the initial conditions are ipo — 0.584 and tq — 5.839. The yellow point denotes the best-fit 

point. 
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FIG. 2: The 68.3%, 95.4% and 99.7% confidence contours with respect to the initial values of (po, ro, using the 
Union 2 dataset. Here we assume the model parameters are ipo = 0.584 and tq — 5.839. The yellow point denotes 

the best-fit point. 

back to we obtain the initial value of po = 1.004. Since the parameter space here is 4-dimensional, we cannot 
show the contour of confidence level in one picture. We have to analyze it in two "cases". First, we fix the initial 
values at their best-fit values, and obtain the corresponding contours of some particular confidence levels of the model 
parameters, as shown in FigH] Second, we fix the model parameter at their best-fit values, and obtain the contours 
of confidence levels of initial conditions, as shown in Figj^l From these two contours, it is easy to find that, to some 
extent, the sensitivity of model parameters and initial conditions has been lowered in torsion cosmology. 

V. FATE OF THE UNIVERSE 

The fate of the universe is an essential issue, which is discussed widely, for almost every cosmological model. In PGT 
cosmology, many works have also been conducted. In [T5l.[T6j. some numerical analyses have been done, which showed 
that H, $ and R have a periodic character at late-time of the evolution for 02 > and 6 > 0, approximately. Some 
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FIG. 3: The late-time evolution of torsion cosmology with the best-fit values of model parameters and initial 

condtions for Union 2 SNela dataset 

follow-up dynamics analysis and statefinder diagnostic done in (2ll . [22j indicate that this character is corresponding 
to an asymptotically stable focus. And the related analytical discussion in [2^ confirmed this conclusion. However, 
the researches mentioned above only presented some qualitative results, where parameters a2 and b and initial values 
are set to certain particular positive values by hand rather than the values constrained via the observational data, 
and thus, quantitative analyses are still needed to be done. Therefore, in order to investigate the evolution of the real 
universe, it is necessary to place the constraints of parameters obtained via the SNela data on this model. 

First we set model parameters and intial conditions to their best- fit values, and solve the evolution equations 
numerically. The solution of Hubble parameter and the trajectory in — $ — i? space are plotted in Fig. |3] and 
Fig. 21 respectively. Then we test some other values of model parameters in the confidence interval of 3ct with fixed 
intial conditions which are still at the best-fit value, and plot the numerical solutions of H{t) in the right column in 
Fig. [5l At last, we fixed the model parameters at their best-fit value, and solve the evolution equations numerically 
with various values of initial conditions in the confidence interval of Scr. Then we plot the results in the left column 
in Fig. [5l From these numerical results, it is easy to find that {H, R) would tend to (0, 0, 0) in the far future, 
which indicate that our universe would expand forever asymptotically to a halt. This description of the fate of the 
universe is in accordance with the earlier analytical work. [2311 
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FIG. 4: The evolution orbit of (H, ^,R) with the best-fit values. The blue point indicates the initial condition 
(1, 0.584, 5.839), whereas the red point (0, 0, 0) denotes the final state of our universe. 




FIG. 5: The late-time evolution of PGT cosmology with various values of model parameters and initial conditions 
lying in the confidence interval of Scr. In the left column, the model parameters are fixed at their best-fit values and 
the initial conditions are varied, whereas in the right column, the initial conditions are fixed at their best-fit values 

and the model parameters are changed. 



VI. SUMMARY AND CONCLUSION 



We studied the cosmology based on PGT with even parity scalar mode dynamical connection, 0+. So the Lagrangian 
we use in this paper takes the form of SNY model Eq.®. We rewrote the motion equations Eqs.(|9|-Eqs. ([TT|) as a 
dimensionless equation w.r.t scale factor a rather than cosmological time. And we obtained the analytical solution 
of past time evolution of this new set of equations, as shown in Eqs. (j23|) -(|25 |) . Then we attempted to investigate 
the goodness of this model, by comparing these theoretical results with the latest observational data. Union 2 SNela 
dataset. Finally, we found the best-fit values of model parameters and initial conditions (02 = 1.336, b = 0.992, 
(fio = 0.584 and rg — 5.839), and that the associated minimal (535.284) is consistent with the ACDM at the la 
level. Furthermore, from the contours of the dynamics analysis conducted in some confidence level FiglT]and[51 it is 
easy to see that, to some extent, the fine-tuning problem has been alleviated in SNY model. Next, we extended our 
investigation to future evolution. We plotted the whole evolution orbit of (H, $, R), from the past to the future, with 
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the best-fit values, in Fig. |4l which gives us a raw picture of the whole evolution. Finally, we tested some other values 
of parameters and initial conditions and found that i/, $, i? all tend to zero in the infinite future, which indicate our 
universe will expand forever asymptotically to a halt. This description of the fate of our universe is consistent with 
the analysis conducted above and the earlier works [2T| - [23{ . Thus, torsion cosmology is a "competitive" model to 
explain the cosmic acceleration, which need not introduce some exotic matter composition. 

In comparison to other models of accelerating universe, torsion cosmology of PGT is new, which still has a great 
number of issues to study. For instance, we could extend the research on SNY model to BHN model, which generalizes 
SNY model to a model with both 0^ and even the coupled term of these two modes. Also, we could investigate the 
effect of 0^ in the very early universe, which might has imprints in CMBR. These issues will considered in the 
upcoming papers. 
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